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ABSTRACT 


The persistence of the fees expansion of the detonation products of a 
high altitude detonation is determined by the sum of the stresses of the en- © 
vironment. Below 200 km the classical dynamic friction with the air couples 
the detonation products 'e6 very much larger mass of air. Above this alti- 
tude more uncertain plasma phenomena of magnetic and electrostatic shocks 
.must be considered to determine this coupling. If a cubic scale height of air 
is raised toa peeateee such that its sound speed is greater ian escape 
velocity, it jets upward, and the ionized fraction is stopped and mixed un- 
stably with the earth's magnetic field. The subsequent expansion along the 
field lines deposits the major fraction of the debris at the opposite conjugate 
point. The B decay in transit on. the various L surfaces reached by. unstable 
mixing tonne the high energy electron injection into trapped orbits. For 
the Starfish event (1.4 megatons at 400 km) picoryooeedicts that the detonation 
products and a cubic scale height of air expands upwards until stopped by the 
earth's field at 1000 km. In addition one expects the deposition of the debris 
in the atmosphere at the southern conjugate point at 250 km.altitude dis- 
persed over 300 km due to magnetic field mixing. Approximately 5% B® injec- 


tion into trapped orbits of L> 1.3 is expected. 
en e-ili- 
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The initial mass-motion of a nuclear explosion is, of course, governed 
by the igrdxcaynarnice of the free expansion of the detonation products; how- 
ever, the subsequent persistence of this motion is determined by the sum of 
the stresses of the environment. Ata high enough altitude the expanding 
detonation products may encounter less than’an equivalent mass of air anda 
free Ricoaeten will take place until some additional stress, as, for instance, 


magnetic field, becomes dominant. 


ATMOSPHERE DOMINATED EXPANSION 


The pertinent question that separates the behavior of the free expansion 
fnotn tae governed by environmental density is: at what altitude will the 
integral mass of air above the explosion that 'collides' with the detonation 
products be greater than the mass of the products ? Since the atmospheric 
density near the earth varies exponentially with height, the supersonic 
momentum- conserving expansion for a snowplow-type shock has a velocity 
distribution 


u(@) = he a) 
| M + i. 4epipje’ dr 
Pease 


“Work done under the auspices of the U. S. Atomic Energy Commission. 
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where 

R = radius of expansion 

= mass of detonation products 

Ug = initial mean expansion velocity 

6 = angle with respect to the vertical 
and 

p(0) = p,e7* 008 8 7 (2) 
where 


h = scale height of the atmosphere 
Pg = atmospheric density (in g/cc) at the detonation point. 


The shape of this shock front at some later time t is 
t 
R(@) = ( u(@) dt. (3) 
0 ; 


If we look al this shape aftcr a time long enough so that R(@ = 0) >h 
‘ 2 cies 
and concern ourselves with cases where the total grams per cm of air is 
large compared to the detonation product mass at a radius large compared to 


h, namely, 
— a ey | 


4nr on. e dr > M, (4) 
0 10) : 


then the radial distribution becomes, Fig. 1 


R(6,t) -u (peace? for 1 - cos @ Pine ; (5) 
0 3 R(0) 
8th Po 


This distribution is based upon the Sasteice that the detonation 
products 'sweep up' the air ahead ina 'snowplow' fashion, i.e., a strong 
shock with the gaiabatic Y - 1. If, on the other hand, the detonation products 
could interpenetrate the ambient air with néplisinic interaction, then very 


obviously the free expansion approximation would apply. 
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DETONATION PRODUCT-AIR INTERACTION 


Three interaction mechanisms between the detonation products and the 
air will be considered: (1) the Coulomb collision dynamic friction (both ion- 
ion and ion-electron), (2) magnetic shocks, and (3) electrostatic shocks. 

Each of these interactions depends upon the initial fractional ionization 
of the air. Two processes for ionizing the air are evident: first, the pulse of 
X rays will give an initial partial ionization, and second, each of the above 
‘mechanisms of interaction will lead to electron heating.and subsequent 
ionization. The energy required for complete ionization represents a small 
eect onal coupling for any of the above three mechanisms since the energy to 
‘be exchanged is very large compared to ionization potentials. In those cases 
where the initial condition for strong interaction requires a modest fractional 
ionization the X-ray pulse will satisfy this condition. 

From Eq. (5), the effective path in which interaction must take place is 
approximately h. The flux of l-kev X rays fees a nominal 1 megaton of 


l-kev X rays at a distance h, assuming small attenuation, is 


b= W,/W,4 n(h)? - 8x10! quanta /em2 (6) 
where 
Wa = 1 megaton 
W_ = 1 kev 
x 


h = scale height = 50 km. 


For a mean absorption coefficient of l=kev X rays in nitrogen and 


¢ 
oxygen of g = 5X ior eae the available energy per atom becomes 


W_ = ooW, = 4 ev. 
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The energy invested per ion pair is usually 32 ev, and considering 
that the Lesombinaion time at the altitudes in question is longer than the 
delay for the arrival of the detonation products At = h/ay ~ 0.015 sec, then 
a fractional ionization of 12% or greater should exist as a minimum initiation 


level for any of the strong interaction mechanisms. 


COLLISIONAL DYNAMIC FRICTION 


ais dynamical friction of a test ion slowing down ina plasma has heen 
treated in detail by Kranzer [Kranzer, 1961 he who confirms the slowing owe 
approximation of Spitzer [Spitzer, 1962]. The dynamic friction on a test ion 
can be separated into two parts — that due to collisions with ions and that due 
to écllisions with electrons. The electron temperature is a controlling fac- 
tor for the electron dynamic friction, because it governs both the ferns 
ization time and the charge state of the ions. (The usual pickup and loss 
dynamic friction for heavy ions passing through neutral matter is smaller 
than the free sivcinon friction for electron temperatures > 5 ev.) | 

The critical region for determining interaction is that region a scale 
height from the detonation where the products arrive essentially noninteracted 
with a velocity Uo and spread in time h/vy seconds. The-maximum heating 
rate then becomes 


a 2 


e 
ate 2No pe Mou jz 
= Z a 
de erg/cc- sec. (7) 


m, aut? /reta,2? gnA 


where No in the number density of the detonation products, M, its atomic 


\ 
mass, Z its charge, and n, the plasma electron density. 


The loss rate of energy from the plasma will be entirely by line emis- 


sion (bremsstrahlung being negligible). An approximation to this loss rate 


C3) 
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can be ‘made according to Post [Post, 1961] and an independent evaluation of 
the Model C Stellarator impurity radiation loss rate [Hoffman et al., 1962] 
can be wade by assuming that each atom which is not stripped below three 
elections has an allowed level 10 ev above the ground eiete with an excitation 


07 16 base 


cross sectiono = 1 This gives rise to a loss rate 


aw /at = non o(kT/2m)!/yy ; | (8) 


5 


= ng. Z(kT/2m)!/2 107! ev/cc: sec. 


An approximation to the ionization state Z can again be made from the 


work of Post as 


pile 


-Z2.% 5< T < 50 ev Bo (9) 


and 


age, 


Letting A, = 50 and 


6 


7 23°). 3 
Nz, = 10° x 6x 10 /Az X 4h 


= 8X 16° ees 
then, equating heating and loss rates, one obtains 


2 m 
an, (Mug jaz | 


Z 2, 1/2fkT\/2 ~15 
1/2 coe, eee Seoicee \eerk 
mi" “(kT)”. /nen To’ enA 


0 


or T= 6x 10° */ng ev, for 1.2 10/9 <n < 1.2.x 1011, (10) 


0 
For Senatice greater than 1.2.x 10'*/ce the approximations (8) for the 
resonant radiation loss sare and (9) for the effective Z of the bomb material 
break down. Since the resonant- state energy is 10 ev, temperatures below 
5 ev give rise to an exponentially reduced paaiuved power, effectively es- 
tablishing a lower limit to the electron temperature of ~5 Sin i addition, 


the pickup and loss cross sections of the detonation products traversing 
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Pa 


neutral matter. at 3 X 108 cm/sec establish.a-minimum effective Z of 
approximately 2. 
At the critical air density, below which there is no interaction,.the 


dynamic friction must be sufficiently large to. slow down the detonation 


products ions ina time t = h/ug: Therefore 


a bie (kT 3/2 
e e 


sec. : (11) 
e nen Zend 


Using .Eq. (9) for n, and Z, and uy = 3X 10° cm/sec and h = 5 
: 
.X 10 cm, then 


2 My m 172 4.3/2 
1.6 X 107° =.—* —€ 


e me ny én A 


where gn A= 15. 


Then 


Ny = 1x iG ce: 


The nuclear Coulomb scattering criterion for detonation producto air 


interaction requires.a still higher density: | 


M, 6x 1071? Ba 2 
n he = =m = 3, 7 a — — = cm. (12) 
U c M c. 2 
air M.. u /2) 
air 0 
For Up = 3.x 108 cm/sec, t.= 6.x 107 )9 cin? 


z* = 50, 2-2 = 600 
aes 
vn, = 1.6% 10!3/ec, 


Thus the electron dynamic friction is considerably larger than the nuclear 


Coulomb scattering. 
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MAGNETIC SHOCK INTERACTION 


The strong interaction by means of a magnetic sheck wave depends 
‘upon the momentum stress of the detonation products compressing.and moving 
the ambient ionized air. If we first neglect the presence of the ionized air 
and ask how the detonation products interact with the field alone, then the 
simple concept of the ionized detonation products blowing a diamagnetic hole 
in the magnetic field is valid. In general, the hole will expand until the work 
done against the field equals the kinetic enesey of the detonation products. 


That is, 


4/3 aR? B*/8r =4X 10°" erg 


R= 10° cm.-for B = 0.3 gauss. (13) 
The thickness of the boundary X between the diamagnetic sphere of expanding 
ionized detonation products and the magnetic field is described by the 'M 


theory' of Rosenbluth [Rosenbluth, 1957]: 


X= ac ie ne2)!/2 (14) 


where m = electron mass, and 


n = number density of electrons associated with detonation 
products at the boundary. 

Since the field can ‘move oe the velocity of light, and the detonation 
Sieducts move very much more slowly, there will be a low danwiky front 
where the momentum.stress equals the field pressure. This is the condition 
for which the above thickness is derived. The physical description of such a 
layer is that the ions penctnate a distance X into the field with negligible 
interaction, but the electrons associated with the ions cannot penetrate the 


field as far because of their smaller Larmor radius. The resulting charge- 


separation electric field decelerates the ions and causes-a transverse drift 
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of the electrons. This electron drift gives rise to a current which is just 
the diamagnetic current bounding the field from the diamagnetic hole. 
For the case of expanding detonation products at a distance h spread 


out in time h/up> the maximum magnetic field interaction distance X becomes; 
2 2 
X = 2500 cm, wheren M uy /2 = BY/8n 


n = 107 for B = 0.5 gauss. 
This then becornes the maximum length that the detonation products can 
penetrate the plasma unless an instability. takes place which causes a break- 
down of the layer. There are two types of instability that can allow the 
plasma to penetrate the field: one is the hydromagnetic Taylor instability, . 
which depends upon the deceleration of heavy fluid by a light fluid, and the 
second is an electron geleciey prawecepace instability in which the electron 
drift stream is unstable to the generation of plasma oscillations. 

The hydromagnetic seta itty, obviously takes place whenever the 
detonation products density is greater than the field-plus-rcoidual.-plasma 
density. As has been shown eucerientaiys [Dickinson et al., 1962], one 
would expect this instability for the conditions of free expansion whenever 
the field stress seesiwanee the detonation products. For the case of the 
earth's field this occurs at R = 1000 kim and so does not affect the initial 
question of coupling to the air. (It will of course, affect the final plasma be- © 
havior that we will discuss ies The éage of coupling toa larger-mass of 
air should therefore be Taylor stable. | 

The sheath instability [Buneman, to be published; Hartman and Sloan, 
to be published; Shkarofsky, to be published] takes place when the electron - 
Larmor radius becomés equal to the Debye length by the generation of plasma 


oscillations which allow the plasma to interpenetrate the field. However, this 
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interpenetration can take place only at a rate given as.a maximum by Bohm 


diffusion D, = ckT /eB. Since this maximum diffusion corresponds toa 


B 
random walk step of a Larmor radius in a Larmor period (ions and electrons), 
- then the effective diffusion velocity must always be small compared to the 
directed stream velocity. Punsdawexts regardless of the instability a strong 
momentum interaction. between the driving detonation products .and the driven 
air plasma will have taken place. The instability merely affects a major un- 
settled problem concerning the shock structure [Gardner et al., 1958; Adlam 
and Allan, 1958; Colgate, 1959; Auer, Hurwitz, and Kolb, 1961; Merawetz, 
1962] namely, the immediate division of internal energy among electrons, 
ions, and ion oscillations. For complete coupling, however, the strength of 
this interaction ieauieestnat either the stress is spread out by instability over 
a distance large compared to'X, or that the ambient magnetic field is com- 


pressed to-a value so that its pressure equals the time rate of change of the 


detonation products momentum. That is, 


ne 
W,/8rh = Bo /8n (15) 
B = 20 gauss. 
max 


This is not too large a compression, buta shock structure of sucha high 


Alfvén Mach number has not been derived on a self-consistent basis. 


NEUTRAL STRESS MAGNETIC SHOCK 


A more conservative picture of the magnetic shock interaction can be 
formed based upon a time-dependent ionization model similar to that formu- 
lated by Gerry, Kantrowitz, and Petschek [Gerry, Kantrowitz, and Petschek, 
1963]. In this model the time-dependent ionization of the neutral fraction oc- 


curring within the moving combined shocked plasma and magnetic field re- 
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sults in the primary momentum stress (Fig. 2). The effective shock thick- 
ness becomes. the ionization mean free path, which may be longer than the 
instability wavelengths discussed previously. As a consequence the shorter 
wavelength instabilities can be neglected, and the ionization mean free path 
in the shocked plasma becomes the critical.parameter that determines 
coupling. The ionization of the neutral air occurs by: 

l. The X-ray pulse (Eq. (6) — the electrons from which may or may 
not have reached equilibrium by the time of arrival of the shock); 

2. The electron ionization; 

3. The charged-detonation-products ions.and shocked-air ions. 

Even assuming equilibrium following the X-ray pulse, this fractional 
ionization is small compared to unity and so the of naey coupling must ie 
determined by 2 and 3. For electron temperatures above 5 ev, the ionization 
rate of air is approximately constant [Mott and Massey. 1949] (o)¥ a3 
x, gage cm * Jae), and so electron ionization should be idendadent of shock 
velovity for air densities where magnetic shocks are important. The lioniza-~ 
tion by partially stripped ions, on the other hand, is sensitively dependent 
upon velocity. Choosing po initial detonation product velocity as the ceitieat 
velocity for coupling, then the equilibrium charge state must be determined 
for the ionization rate. Knipp and Teller [Knipp and Beiter: cae have | 
satiated the effective charge of a ae zZ atom slowing down ina Heueral gas 


and have shown that capture and loss are in equilibrium when Ug = Vy Ties. 


where Ye is the bound orbital electron velocity. Choosing Up = 3X 10° cm/sec, 
then the ionization energy corresponding to ve 4.5% 108 cm/sec is 80 ev or 
Z ost = 3.5. 


The cross section for nitrogen ionization by an ion of charge Z off and 


velocity Uo is 
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«6 10°° (Z, 2S gig. ea 
. The detonation product density at a distance h is 


Ny = Mx 6x oo, 4/3nh> For -M = 10° gram and A = 20, 


ND 


3x {0 Jee: . . (16) 


Therefore the mean free path of a neutral air atom in the flux of detonation 

products becomes | 
d = 1/No 356 10° ean 

The electron density such that electron ionization will have contributed 

equally to the ionization of the air within the time of the passage of the detona- 

tion products is: 


9 


ns = Up [ovr = 2 X10’ electrons/cc. _ | an) 


e 

In summary, a strong interaction between detonation products and the 
combined air and magnetic field is predicted for an ionization-type hydro- 
magnetic shock where the shock thickness due to either ions or electron 
ionization is large‘ and comparable to the scale height of the atmosphere. 
During the expansion phase such a large shock thickness should be stable un- 
til the magnetic field stress reverses the expansion. 

Since the ionization magnetic shock by debris ions becomes the dominant 
shock mechanism above 350 km, it is pertinent to demonstrate that the high 
aeuhts charge state of 3.5 may have had a chance to reach and maintain 
equities before penis expanded to a scale height in radius. It is evident 
that a free expansion in a vacuum will permit the detonation. products to cool 

sy adiabatic expansion an radiation leading to recombination and a lowering 
of charge state. ever: because of collisions with ait molecules the 
detonation products will be heated, but this heating is small compared to 


‘blackbody radiation (even at 1 ev). Consequently, it is only: after the detonation 
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products have saenaed to fhe point where they are optically thin, yet an air- 
atom collision range thick, that the temperature will rise and the anes 
state become high. Naturally, if the detonation panies are a collision 
range thick, there is more than ample time for charge equilibrium to be 
reached. | 
The range from.Eq. (11) becomes 
NpR =5*X ion! atoms /em* _ 


where 


Ny = detonation product atom density 


R = radius of detonation products. 
Solving for R analogously to Eq. (16) gives 
R= 1km, Ny =5-X 1014 fee. 
Using Eq. (8) for Ane oaGuank line emission power loss (an upper limit 
because of the available states and self-dbacrpion)and assuming 
n= 3.5 Ny gives a power loss rate 
12 


2 
Pos 10 ergs/cm ‘sec. 


- The air density Pp. such that the collisional heating equals the power radiated 
Prad ~ Pa ny 2 | 
becomes | 
subir Pees 
At lower air densities the balance will occur at a larger detonation product 
radius since the heating rate decreases as 1/Ny for RS air sheen range; 
whereas the radiation decreases as 1/N Consequently, over the air- 


density range of interest ( igo? SP io." 


= the detonation-product charge 
state should be maintained at an equilibrium value of Z = 3.5 until slowed 


down, - 
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ELECTROSTATIC SHOCK INTERACTION 


i 


The electrostatic shock pereer upon the growth Sf lange acapinvade 
plasma oscillations and /or ion sound waves due to the two-stream interaction. 
Buneman [Buneman, 1959] has shown that the electron-ion two-stream in- 
stability grows .to saturation in roughly 30 plasma periods. Saturation occurs 
when the electron stream kinetic energy is randomized by noniiness effects 
to a quasi-thermal distribution. When the resulting ieetten thermal velocity 
(Vern) equals the relative stream uetocies linear-unstable growth stops. The 
ion-wave instability as shown in detail by Stringer [Stringer, 1961] grows -pro= 
vided the ion temperature is small compared to the electron temperature 
and provided the electron thermal energy is less than the ion stream energy: 
On the other hand, the maximum ion wave amplitude (measured in neve) is 
limited in linear theory to the electron temperature. The possibility of 
further ‘nonlinear instability growth has been investigated numerically by 
‘Hartman and Colgate [Hartman and Colgate, 1962], and although large 
scnelitade growth has been confirmed for electron ion mass ratios up-to 100, 
sorne doubt exists for-mass ratios as gheat as 10° for the present case. 

ew: Reventine the ‘more sececweee result of the saturation of the electron 


instability at Voth = Yo? this gives a two-stream dynamic friction that will 


th. 
maintain the electron temperatures at m re /2 = 50 ev. This magnitude of 
the friction is then limited by the plasma radiation loss rate at 50 ev. 


Equating this loss rate to the detonation products energy input rate determines 


a minimum density for interaction. From Eqs. (8) and (9), 


z ii plz (kT /2m)!/2 o(hv)k = Wu,/8nh* . (18): 
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where 


10 ev, (eT /2m) 1/2 =U, = 3x 108 cm/sec, T = 50 ev, and 


= 10. 16 eine 


a 
z 
il} 


Then 


n, = 1.4*X iO fee: 


This is a stronger interaction than the electron Coulomb dynamic friction, 


but certainly weaker than the magnetic shock. 


ENERGY CONSERVING INTERACTION AND EXPANSION 


The magnetic shock ensures strong interaction essentially independent 
of density, but the limit of free expansion behavior is determined by the radial 
distribution of Eq. (5). Choosing M = 10° grams, and h = 50 km, then Pain 
becomes}! 


=M/8rh? = 3x 107? gee . (19) 


Pmin 


or No * 1.25 X fo fees 


Above this altitude free expansion takes place. Below this altitude the com- 
bination of detonation products and heated air will expand upward at a velocity 
given by Eq. (1) provided the resonant-line radiation epeting of the air behind 
the shock justifies the approximation of a snowplow shock, that is, that the 
internal energy is radiated away as fast as it forms. However, for densities 


less than the limits of Eq. (10), that is, io*® <n, < 10' /ee, the resonant 


0 
radiation loss very rapidly becomes negligible, because air atoms are stripped 
of electrons below the lithium 3-electron configuration, and only levels of 
large energy gap and small cross section remain. The radiation: power loss 


rate from air decreases by two orders of magnitude as the temperature in- 


creases from 20 to 100 ev [Post, 1961]. Consequently, in this intermediate 
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density region the subsequent expansion of the high-temperature air plasma 
‘controls the hydrodynamics more than the directed Hemet an of the detonation 
products expansion. 

In order for this 'bubble' of high temperature plasma.to expand up- 
wards through the remaining atmosphere, it must comprise a volume of 
roughly-a cubic scale height, otherwise the mas oy above will be greater 
than that in the 'bubble', thereby limiting upward motion. The expansion 
velocity upward is then sound speed within the bubbles Assuming a specific 
heat nee electron of 3 to include the er of eer and the radiation during 


ionization and a volume frorn Eq. (5), the expansion velocity A. becomes 


3 2 /> 
24nph™ (u.,) Jes Ww, 


ar 1/2 
or Uo = (w, /l2np | cm/sec. (20) 


The bubble will cool by both radiation and expansion, but in general 
once the temperature has run away above the resonant line-emission range, 
the adiabatic cooling will dominate. As pointed out earlier, the subsequent 
Spansion will be limited by the stress of the magnetic field at roughly 1000 


km for a megaton. Assuming a Y = 5/3 adiabat, the final temperature T; be- 


comes 

Bay ie WR (21) 

e f t 
whece n =8.2n. for T, > 50 ev. 
e 0 i . 
Then 
\ : 2 
Pe w,/25 ny k Rv h. (22) 


For R= 10° km, h = 50 km, Ww, =] megaton, and 


— q 10 a 
T, = 18x 10 {a9 ev. 
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_ When the bubble material is stopped by the magnetic field, one would expect 
" some heating ae the electrons from the inelastic energy of the ions. The 
‘fractional cooling to be expected within the time of interaction between the 
bubble material and magnetic field (assuming a lies within the resonance 


radiation range) by adiabatic cooling is by Eq. (8): 


ae (h/R)? n, (2K Tn)? (hv yk R/u.y, 


a 3 8 
=n, T/us, 1.6X10-. (23) 


Z 


This fraction is small (~ 10°“) for the largest No (~ 10 Jee: and so the 


radiative cooling of the plasma after expansion is negligible. So also is the 
heating for T> T « T is that electron temperature such that the 
- max max 


electrons can make one thermal exchange time with themselves within the 


expansion time. 
l/n. ove R/u 
/ e / cx 


nen ti/2 n/a), o = 6x 10713 /7? om (24) 


0 


< 
i 


= 5.5 10! gil cm/sec 


T=0.3n Ju ev. 
0° “ex 


When the adiabatic temperature of Eq. (22) exceeds the heating temperature 
given by Eq. (24), then no further change should be expected, At Ny = 10°" /ce 


the electron temperature should stabilize at roughly 10 ev from both effects. 


MAGNETIC FIELD LIMITED EXPANSION AND FINAL DEPOSITION 


When the plasma of air and detonation products is stopped by the stress 


of the magnetic field at a distance R, we expect a Rayleigh-Taylor unstable 
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mixing between the two because the density of the incident plasma is greater 
than :the ambient density (Fig. 3). The linear growth of a perturbation of 


amplitude Ay .is [Dickinson et al., 1962]: 
A= Ay exp [g./2m) [2 cal (25) 
where g is the acceleration and \ the wavelength. For a linear deceleration 


taking place in a distance R, the growth becomes: 


A= Ay exp IRAr 2, (26) 


Only the soniiveds phase of instability growth (where.A> \/2) will irreversibly 
mix the plasma with the field. Consequently, if we wish to know the irre- 
versible mixing of the plasma onto magnetic flux surfaces other than the one 
defined by the explosion, we must ask what is.the largest wavelength that re- 
quires the deweut generations of growth to reach an amplitude greater than 
/2. Then, because of the faster growth rate of smaller wavelengths, all 
smalle} wavelengths should have reached the Aontineas limit too. 

For a rising bubble or jet of plasma of form approximated by Eq. (5) 
the wavelength \.= Rf has an amplitude of approximately \/2 so that during 
deceleration by the field by Eq. (26) one generation further growth can take 
place within the nonlinear limit. Therefore, the plasma should penetrate 
the field with a mixing length R/r. 

The subsequent motion of the plasma will be an expansion along the 
lines of fsxee, since in the parallel direction less stress is exerted by the 
field. The principle direction of this expansion will be away from the ore 
aches of detonation because of the following stresses, all of which tend is 
direct the plasma from a northern hemisphere Suisse towards a southern 


hemisphere deposition: 
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1. The verticle motion of initial expansion Eq. (5) has a component 
of momentum along the field lines towards the far conjugate point. 

2. The ambient field initially diverges in a direction towards the far 
conjugate point. 

3. The distorted field diverges more strongly because of the 'line 
tying’ in the aetialiy ionized atmosphere at the near conjugate 
point. 

For these reasons the plasma will travel along the lines of force until they 
intersect the earth's atmosphere at the far conjugate point with a distribution 
determined by the magnetic flux surfaces intersection with the earth. Since 
we have already shown that due to the low ages of the plasma at the time 
of mixing the interaction should be energy~-conserving, the velocity along the 


lines will be sound speed, namely, Yo 


DEPOSITION IN THE ATMOSPHERE AT THE CONJUGATE POINT 


On intersection with the atmosphere at the conjugate point, the plasma 
will have a width (longitude) as well as latitude determined by both the mixing 
length R/n and the diamagnetic pressure. Since the plasma expands intn a 
tube roughly 6000 km long, the radius for a pressure the same as the bubble 
of R = 1000 km is also approximately 300 km, The altitude of deposition will 
be determined by the electron dynamic friction and the nuclear Coulomb 
seattering. The impact waiodiyy Vox determines which of these two ay aaa 
frictions will dominate. If the atmosphere is not heated by impact toa 
temperature high enough so that the free electron density is comparable to the 
nuclear density, then the electron dynamic friction becomes small for an 
impact Seleeiy ‘such that the relative kinetic energy of an electron is less 


than the binding energy in the atom, namely: 
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aes /2> Tew 
ex = 
or 


8 . 
ed 1.5 10° cm/sec. (27) 


Below this velocity the bound-free electron dynamic friction becomes 
small and so-the heating of the atmosphere is small until the nuclear scattering 
becomes large. By then the incident plasma is stopped and the question of 
heating the atmosphere to ionisation becomes unimportant. Only in.the 
velocity range as 1.5.X 10° does the question of heating become important. 
The energy available is reduced to w,/3 due to radiation loss during’ ioniza- 
tion and recombination, and this energy is deposited over an area (R/3)°. 
Assuming twice the ionization potential in electron volts is required to ionize 
each atom because of radiation loss, the maximum density Nass that could 


be ionized is: 


2 
hN vay 2 Ke W,/(R/3) . ' (28) 
for 
E. = 7 ev. 
1 ‘ 
N = 1.2.x 0 ew. 
max : 


This is too low a density to stop the plasma by electron free-free 


dynamic friction by Eq. (11), so the dominant stopping will be by bound-free | 


8 8 


for og? 1.5.X 10° cm/sec and by nuclear-Coulomb for Uo <1.5.X 10°. 


The nuclear-Coulomb dynamic friction is given by Eq. (12) on the basis that 
both nuclei of a collision are unscreened by bound electrons. This is the 


1/2 


case if the impact parameter (c/r) is small compared to the orbit of the 


k-shell: electron. 
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‘The minimum incident energy for which Eq. (12) is valid is determined 


oe 
Aen _ An a 2+ on A , (29) 
(Ry P2777) 2 (M_. a? /2)* | 
” air “ex 
where Ry = Rydberg constant = 13 ev 
and j gnAz= lo 
giving 
Mu /2 = 13 kev. 
exp 


At lower incident energies the term that sums the multiple scattering 
(¢gnA) approaches unity and the cross section approaches the classical 
Rutherford cross section at M us /2=4kev. Therefore an approximately 
constant stopping cross section of o = 1.5X foo eas can be given for the 
energy range 4 < Mu” /2 < 13 kev. This corresponds to deposition at an 
altitude of 250 km. | 

Since for the very much higher energy region the bound-free dynamic 
friction should stop the incident plasma at the same altitude (Eq. (11); 
Yes = uy /2), we can feel reasonably confident of a deposition height of 250 


km. The distribution according to magnetic flux surface should therefore be 


compared to the earth's field at this altitude. 


MAGNETIC FIELD GRADIENT DRIFT 


In addition to the spreading in flux surface due to the Taylor instability ) 
_mixing, there will also be a drift in longitude due to finite Larmor orbit drift 
-in the gradient of the magnetic field. This depends upon the following 


parameters and assumptions: 
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1. The vegion of most effective field gradient may be the small region 
of high gradient in the initial field distortion or the dipole gradient effective 
during the transit to the conjugate point. 

2. The transverse electric field established by charge separation may, 
or may not dominate the VB drift. 

3. The transverse ion velocity and equilibrium ene state effect both 
the EX B and VB drift. 

Unfortunately, both regions of field — the distortion loop and dipole 
gradient — are Sevive in producing drift withthe greatest uncertainty con- 


cerning the detailed loop shape and gradients. In addition the relative im- _ 
portance of the electron stress versus ion stress in the formation of the vB 


of the loop is dependent upon a detailed knowledge of the nonlinear instability 
mixing behavior. If the interaction boundary is a stable 'M' layer as al- 
ready discussed, then the ion stress is balanced by an electric field from 
charge separation and only the electrons drift transverse .to-the magnetic ; 
field. Onthe other hand, if the layer is unstable, the stress will be divided | 
according to the respective temperatures — which in turn are dependent upon — 
the form of instability. 

The drift in the dipole field and in the distortion loop will be calculated 


on the a priori assumptions: 


I 
\ 


(1) The electron temperature is determined by the linear processes 


only and will be taken as 10 ev from Eq. (24) for the limited initial dens ity 


10 


region 10? < ny < 10° 


(2) The ion transverse velocity in the earth's ficld is the expansion 


velocity u__. 
ex 
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For these assumptions the drift velocity YD becomes 


Vv =u 


= neta VB /B 


r, = ion Larmor radius =Mu_ /Z B. (30) 
ex’ “e 


L 
The charge state Z is determined by the recombination ee in the 

expanded plasma at the density Ny (h/R)? and T. ~ 10 ev. This rate (for 

3-electron configurations) is given by Post [Post, 1961 ]-as 

agx 107124 


nov.=n oF sey) (31) 
3/2. 3 
Te, , . 


where k,(y) will be close to unity for kT) << E., the ionization potential of. 
the state Z off 


For a recombination time of 1 sec, 


3,1/4 
= 1.7™x 
Zoe = 1-7 X 10 /n 
for T: = 10 ew 
e 
giving Zaee * Lod. 


This means that for the anni range considered, radiative recombination 

is much too slow a process to occur within the available transit time of the 
plasma. Consequently, the charge state will be determined by the eons 
charge stripped state during the initial detonation products-air interaction. 
This in turn depends upon the runaway electron temperature in .this region - 
and therefore upon the form of the interaction. ao Eqs. (9) and (10) we 
will assume the air has been completely stripped to a mean charge of 7. The 
subsequent charge exchange with neutral air atoms in traversing a distance 
of the.earth's radius at an ambient density of 10° /ee is small. From Mott 


and Massey [Mott and Massey, 1949 | 
a 2 
o = 10 16 Co /2.5 x 10°) where Me 2.5™& 10° cm/sec. (32) 
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The fraction charge exchange to Z = 6 is less.than 


it ie eb 1074 


(L = earth's radius = 6X 108 cm) 
and so is negligible. 
The Larmor radius Men becomes 


r, =5.3X 1074 /B ee car eee ee field. 
L ; ex ex 


To calculate the drift in the earth's field the dipole field gradient is 
VB/B = 3/1 
so that the drift velocity becomes: 


2 -12 
Vp = ee 3.5.% 10 cm/sec, 


‘and the drift D becomes 


7 -3 . 
D = Lv )/4.,. = 2x10" ul cm. (33) 


Some fraction of the detonation products as a precursor may not have been 
subjected to the runaway temperature region of the mean interaction and so 
may have a much lower state of ionization. Choosing the extreme of charge 
one for a fission product the drift becomes 

Dea oem, p 3 (34) 


ex 
This distance is smaller than the mixing length spread provided ee < Uo 
and so should not be observable.: | 

The drift in the formation of the distortion loop may be considerably 

larger. Again assuming the stress oh the field is created by the ions, then 
the logarithmic gradient of the field should be roughly the ion Larmor radius 
_ for a length of time corresponding to the intermixing of plasma ree field. 
‘ Since the drift velocity becomes roughly half the incident velocity Wey? and 
the period of drift is R/t Wox? the total drift becomes. 


D x R/er. . . . (35) 
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BETA INJECTION FROM RADIOACTIVE DEBRIS 


The injection of energetic electrons from the B decay of the radio- 
active debris should occur daviay thé transit of the debris to the conjugate 
point. Only the decay after the unstable mixing of ae plasma and field should 
permit injection onto flux surfaces other than the one intersecting. the detona- 
tion point. Since the time to reach the unstable eee oe condition R/u is in 
general shorter (< 1 sec) than the major fraction of the B decay (~3 sec), and 
the flight time to the conjugate point longer than 5 seconda,: then most of the 
decay 6's will be injected into flux surfaces above tne atmosphere. However, . 
eng those flux surfaces that do not intersect the earth's atmosphere for a 
reasonable mirror ratio will show long ® lifetime. This minimum L surface 
[McUwain, 1961] is approximately L = 1.3 so that only that fraction of the 
debris reaching this L surface and above:by unstable mixing should contribute 


to the long-lifetime 8 injection. 


NEUTRAL AND SINGLE-CHARGED DETONATION PRODUCT BEHAVIOR 


The considerations so far have-been for an average behavior and only 
the dominant interactions have been considered. It is evident, however, that 
the charge state will Austad and the effects of a statistically small fraction 
| may be important. In particular, the singly-charged high energy tail of the 
detonation dudduets may drive the magnetic field and partially (X-ray) ionized 
air away from the remaining neutral fraction bevore ionization is complete, 
and secondly, a small fraction of the debris neutralized by nideee exchange 
will expand across the field lines to give B-decay injection at higher L values 


than would be reached by the Taylor instability mixing. 
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BETA INJECTION BY NEUTRAL DEBRIS 


The fatten of the debris that escapes as neutral atoms should be 
generated predominantly by the slower ‘moving depuis that expands upward 
after intersecting with the air. The higher velocity precursor debris has a 
smaller puShebiiy for escaping as neutrals. The ee section for pickup. 


.from neutral air is [Mott and Massey, 1949] 


8 5 
sete 10 2 8 
19 = 10 (2 cm for err 10° cm/sec (36) 
ex 
= 10. eas aa for u < 108. : 
ex — 


On the other hand, the air through which the debris is rising and expanding 
will be inived by the debris. Only that density region which is not fully 
ionized by the expanding debris will contribute to the charge neutralization. . 
Assuming the ae temperature during this phase of expansion is greater 
than 5 ev, the ionization rate will be temperature-independent and the ex- 


panded density determines when the further ionization of ambient neutrals 


is unlikely. This density is determined by 
n o,v=zu_/h, (37) 


e i “ex 
giving 
n =5u 
e ex 


The expansion of the plasma at initial density 6 No to the density n, 


takes place ina distance r into a solid angle of 1 radian such that: 
: fs | 
r=h(3 X6n)/5 u_) (38) 


and the ambient neutral density becomes: 


pan 


No = Ny exp +“ rol ted | ; (39) 
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‘The neutral fection becomes 
f= No Toh ~ (40) 


7 -9 1/3 is 8 
=5%X 10 Ny ©XP [8 Ny /a,.,.) | tox Wao 10°. 


This is a major fraction only for ny ® 108, or an altitude greater than 550 ‘Am. 
The £§ injection from this deubvanived debris will occur according to 
the usual t_ ie power law for fission product decay. It should be noted that 
for fairly low altitude detonations where Vos is very small.(Eq. (20)) — of the 
order 10° cm/sec — neutralization will take place by recombination within 
the expanding bubble. This debris will be rising slowly so that both p decay 
and sunlight will cause beteatoaticn, The decay spectrum will change within 
this time to give a different injection spectrum as a function of L value. ‘In 
addition, a significant change in.the trapped £ spectrum should be abserved 
for flux surfaces reached in. times of the order of an hour, and the f velocity 
distribution function at the equator will depend upon the debris reaching | 


their corresponding geopotential. 


AURORAL EFFECTS 


.The deposition of the detonation products and air plasma at the far 
conjugate Senet heat up the air in which it is deposited, causing the air 
to expand upward and intersecting some of the trapped B-decay edtabna. 
The optical pattern will correspond to the illumination of the flux'surfaces 
on which high-energy electrons are trapped and hence to the gradient field 
drift separated charge state of the B-decaying debris. 

As previously discussed, the heating takes place at 250 km, depositing 
Wo electron volts per atom during the transit time of the plasma. This 


energy. is: 


’ 


D 
for R/3 = 300 km, n 
k 


then Wy 


= 40 ev/atom. 
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Wp = W,/3N, = 4x 1077/3 WR/3)° ak Sa (41) 


2.x 19/9 and 


1.6.X tor 7 ergs /ev, 


The transit time is so long, however, (up to tens of seconds) that this energy 


will be almost entirely radiated away as it is deposited. The residual 


temperature will therefore be determined by the radiation properties of the 


air — which becomes transparent at roughly 1/2-ev temperature. The 


_vertical expansion velocity corresponding to twice sound speed at this.tempera- 


ture will be approximately 3 X 10° cm/sec, permitting the air to reach 500 


to 1000 km.altitude in several minutes. The optical radiation corresponds | 


to the excitation of the air with approximately 10% fluorescent efficiency for 


an energy deposition rate corresponding to the usual minimum ionization 


‘rate of 2 Mev /g/em”. 


The 6-decay electron density for uniform phase-space injection into a 


magnetic tube of force 300 km wide and 6000 km long for an injection time 


longer than 5 seconds becomes: 


n.= W/E. Vol = 2.5X 10° electron/cc - (42) 
“eo 6 ie : : 
where 
We = B-decay energy = 0.05 Wy 
E = mean B energy = 1 Mev, and 


Vol = Volume of flux tube. 


The radiated energy loss rate then becomes 


R 


6 


Ey = engp X 2X 10° = ld" “eee ie sabe | (43) 
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where c = electron velocity 
ts = aenety after -10-fold expansion 
| _=5®& io g/cc 
= fuceseeent efficiency z 10%. 
"The total optical brightness B. fora peneGawent thickness of 300 km 


i 


‘becomes 


Bae x 10! E_ = 30 exg/em’* see. (44) 


The heating rate of the air is ten-fold bigger than Eq. (43) and corresponds 

' tol ev _— atom during the 100 seconds of expansion. This is enough to 
enhance the: eepansion so that some air would be expected to ionize and expand > 
Baek along the lines of force from both conjugate regions.to meet at the 
equator. - The energy loss rate per electron is 2000 ev/sec, but during the 
available time before the gradient field drift carries the electrons eastward 
the fractional energy loss is small. . The energy loss aceurs for a time: 


Ts (R/3)/V, = 30 sec where 


Vp 7 10° cm/sec from Eq. (27) for a 1-Mev electron 


so-that only the low energy electrons with larger energy loss rate and slower 


drift will be lost in the rising air mass. 
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Fiz. 1. The time sequence of the expansion of detonation products ina stratified atmosphere is 
represented by the approximate snowplow shock contour. The air and detonation products in this case 
are considered neutral and only subsequent fission decay gives debris interacting with magnetic field. 
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Fig. 2. “Zonization magnetic shock (moving frame). 
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- LEGAL NOTICE 


This report was prepared as an account of Government sponsored work. 
Neither the United States, nor the Commission, nor any person acting on 
behalf of the Commission: 


A. Makes any warranty or representation, expressed or implied, with 
respect to the accuracy, completeness, or usefulness of the information con- 
tained in this report, or that the use of any information, apparatus, method, 
or process disclosed in this report may not infringe privately owned rights; or 


B. Assumes any liabilities with respect to the use of, or for damages 


resulting from the use of any information, apparatus, method or process dis- 
closed in this report. 


As used in the above, "person acting on behalf of the Commission " 


includes any employee or contractor of the commission, or employee of such 
contractor, to the extent that such employee or contractor of the Commission, 
or employee of such contractor prepares, disseminates, or provides access 
to, any information pursuant to his employment or contract with the Commis- 
sion, or his employment with such contractor. 


